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Abstract. High-resolution observations show the fine structure of the global equilibrium magnetic 
field configuration in solar atmosphere to be essentially different from that assumed in the tra- 
ditional "potential + force-free" field scenarios. The interacting large-scale structures of fine field 
elements are separated by numerous non-force-free elements (tangential discontinuities) which are 
neglected in the traditional field picture. An incorporation of these elements into the model implies 
a dynamical rather than statical character of equilibrium of the field configuration. A transition of 
the system into flaring can be triggered by the ballooning mode of flute instability of prominences 
or/and coronal condensations. Tearing-mode and MHD instabilities as well as the effects of over- 
heating of the turbulent current sheet prevent the field from stationary reconnection as it is adopted 
in the traditional scenario. We speculate around the assumption that the energy release in active re- 
gions is governed by the same scenario as dynamical current percolation through a random network 
of resistors in which the saltatory resistance is controlled by local current. 
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INTRODUCTION 

High resolution observations of solar atmosphere with space telescopes (e.g. SOHO, 
SDO, TRACE, HINODE, etc.) revealed that some parameters of the solar flares are 
inconsistent with the predictions of currently used theoretical models. This reopens the 
following questions: 

• What is the equilibrium configuration of the magnetic field and currents in the 
active regions before the flaring events? Is this a static or dynamical equilibrium 
and in the latter case, what is the topology of the phase trajectories attractor in the 
multi-dimension phase space of the system? 

• What could be the trigger of flares? Is the equilibrium system stable with respect to 
small-amplitude external perturbations? 

• How does the field and current structures evolve during the flare process and what 
is the mechanism of energy release? What is the post-flare configuration of the 
magnetic field? 



Possible answers to these questions are discussed in this paper. 

PRE-FLARE FIELD CONFIGURATION 

Pre-flare equilibrium field configuration is usually treated in terms of a so called 
"potential+force-free paradigm" [18, 27]. The field configuration within this approach 
is considered as a superposition H = Ho + Hff, where Ho is the potential component 
(AHo = 0) and Hff is the force-free component of the magnetic field, which satisfies the 
condition rotHff = aH ff . Parameter a is accounting for the field helicity and in the gen- 
eral case can be inhomogeneous across the force tube of the magnetic field. The potential 
field can be unambiguously calculated if the boundary conditions (three components of 
the field at the photosphere) are specified. The task of determining the force-free com- 
ponent is more complicated. Only a limited number of analytical solutions have been 
derived for the most simple field configurations and under the assumption on an infi- 
nite half-space. Progress has been achieved in 2D numerical simulations. The boundary 
conditions in these calculations have been taken from observations and only a few sim- 
plifications about the structure and helicity of the magnetic field in the active region have 
been invoked (see e.g. [25], and references therein). A relatively good reproduction of 
the observed low-resolution maps of the magnetic field has been achieved for certain 
values of the free parameters of the model. Incorporation of more complicated boundary 
conditions (e.g. two or more dipoles or/and a dipole field in an external magnetic field) 
into the model has led, however, to singularities in which an unambiguous determination 
of the magnetic field is impossible. These singularities in 2D calculation with non-static 
boundary conditions have been associated with regions in which the magnetic energy is 
accumulated and released during flares. 

The above conclusion has been seriously challenged in studies of 3D configurations 
in which the tangential component of the magnetic field is taken into account. The 
singular points in this case do not appear and the field configuration can be explained in 
terms of a magnetic field with shear, which is similar to the field configuration realized 
in TOKAMAK [7]. A formation of a global static force-free magnetic configuration 
resulting from the interacting local force-free field structures is in this case impossible. 
Instead, the interaction between the local force-free field regions leads to a formation 
of tangential discontinuities [16] and re-construction of the global field configuration 
of a dynamically equilibrium state associated with the attractor in the phase space 
[jq,Xi]. This clarifies why the magnetic energy is accumulated in singular points in 2D 
simulations. This is a consequence of a basically incorrect assumption about the zero 
value of the magnetic field in the vicinity of these points and, correspondingly, the zero 
value of the magnetic tension (the energy of waves propagating in the vicinity of these 
points accumulates due to an infinitely large value of the refraction coefficient). 

Analysis of high-resolution images (see e.g. the images of the sun taken by TRACE at 
http://soi.stanford.edu/results/SolPhys200/Schrijver/TRACEpodarchive.html) suggests 
that the solar atmosphere is threaded by a large number of narrow magnetic flux tubes 
connecting the chromosphere and the corona. The radius of the tube cross-section 
ranges from a few tens to a few hundred kilometers and remains almost constant along 
the tube from the bottom to the top. As the magnetic flux in the tubes is conserved, 



the field strength at the top of the tubes in the corona can be as high as hundreds (or 
even thousands) of Gauss. This indicates that the magnetic energy transferred from 
the photosphere to the corona is concentrated predominantly in the numerous narrow 
magnetic tubes which represent a global magnetic structure in dynamical equilibrium 
[17]. This significantly differs from the picture expected within the 'potential+force-free 
paradigm' (see e.g. [26], [15]). 

Observations of the fine structure of the magnetic field in the solar photosphere [29] 
show some evidence for the fractal structure of the magnetic field in active regions. 
This favors scenarios in which the fine structure is an intrinsic property of the magnetic 
field emerging from the convective region rather than a product of instabilities or/and 
turbulent motions in the solar atmosphere. Resulting pre-flare magnetic structure is an 
ensemble of numerous strongly interacting thin magnetic threads with avalanche-like 
percolation of free magnetic energy through this complex network [31, 1, 14]. 

FLARE TRIGGERS AND ENERGY RELEASE PROCESSES 

Turbulent current sheet paradigm 

Solar flares are widely associated with anomalous dissipation of magnetic energy 
in the turbulent current sheet. Turbulization of the material in a region with strong 
current can be expected if the velocity of electrons exceeds the phase velocity of the 
plasma waves M curr = j/ne > V p h, where j is the current density, and n and e are 
the number density and electric charge of the electrons. Otherwise, the waves will 
be suppressed by the Landau damping. The phase velocity of the plasma waves in 
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the case of the ion-acoustic turbulence is V p h = (T & /m\) 1 , where T s is the electron 
temperature and m x the ion mass. Hence, a formation of a turbulent current sheet in the 
solar atmosphere (H ~ 100 G, n ~ 10 8 cm - 3 ,andr~10 6 K) could occur if the thickness 
of the current layer were 1-10 km, i.e. a factor of 10 3 — 10 4 smaller than the size of the 
active region itself . For the current sheet to remain active during the flare in ~ 10 10 cm~ 3 , 
T ~ 10 7 - 5 K), its thickness should not exceed 1-10 cm. Under these conditions, however, 
the current sheet would be disintegrated by the tearing instability in a very short time. 
A question about the stability of the current sheet on the timescale of flaring events, 
therefore, remains open so far. 



Trigger mechanisms 

The observed rapid transition of the system from its equilibrium state to a flare in- 
dicates that the magnetic configuration in which the energy is stored before the flare is 
metastable. In other words, the system before the flare is in a pre-catastrophe equilibrium 
in which a slight external perturbation, caused by the instability of a surrounding struc- 
ture, may act as a trigger of flaring. Among the magnetic configurations which meets this 
criteria are the prominences and the coronal condensations. As shown in [19, 20], these 
configurations are interchange unstable. In particular, a transition of these configurations 



into a flare can occur as the result of a flute instability. Observations of pre-flaring states 
have revealed that: (i) Almost all flares are preceded by emerging of the magnetic flux 
which supplies the magnetic energy into the active region [4, 5, 6, 13, 24], (ii) Some 
flares were observed to be preceded by oscillations of prominences in the vertical di- 
rection, expected at the first stage of the flute instability development [12, 19], (iii) The 
strongest flares have been observed to be preceded by raising oscillations in coronal con- 
densations [9, 23, 22], and (iv) The flares develop as the phase wave propagates through 
the magnetic arches, starting from a local region of initial energy release. The magnetic 
energy in each of the arches is converted into plasma heating and particle acceleration up 
to relativistic energies. The energy of the particles propagating towards the solar surface 
is then converted into heating of chromospheric gas which fills the magnetic tubes and 
flows up to the corona, forming the coronal mass ejections. Thus, imbalance of unstable 
magnetic configurations by an additional magnetic flux emerging from the convective 
region may lead to a development of interchange instabilities of the prominences and 
coronal condensations and trigger the phase transition of the global magnetic structure. 

ENERGY RELEASE IN THE CURRENT SHEET 

It is presently well established that flares start with the formation of turbulent current 
sheet [8]. However, the life-time of the current sheet is limited due to the following 
reasons. 

Current sheet instabilities. The characteristic time of tearing instability in the current 
sheets of thickness d, can be evaluated as T tg ~ ^A T d ^> wnere t a = d/V& is the Alfven 
time and V A the Alfven velocity. Td = Anod 2 /c 2 is the diffusion time and a is the plasma 
conductivity in the current sheet. Finally, the value of s depends on the type of instability 
and lies in the range 0.25-0.5. The tearing instability of the current sheet is followed by 
pinch-type MHD instabilities (sausage and kink modes) on a time scale T p j nc h ~ Ta, 
which leads to the disintegration of the current sheet. 

Overheating of the current sheet. A flare starts as turbulization of the plasma in 
the current sheet occurs. The onset condition for turbulization reads j > j CI , where 
jet = neu CI is the critical current density and w cr the critical velocity, which is equal 
to V Te = (kT/m e ) l/2 for T e = 7- = T, or c Si = (kT e / mi ) 1/2 for T e > 7- [20]. Here n is the 
number density and T e and T\ are the electron and ion temperatures in the current sheet. 
As the turbulization starts, the temperature of plasma in the current sheet increases on a 
time scale of ~ 1 /oooi, which under the conditions of interest is about a few microseconds 
(here (% is the ion plasma frequency). If the cooling of the plasma in the current sheet is 
ineffective, the gas temperature reaches a critical value at which the onset condition for 
plasma turbulization is no longer satisfied and plasma waves are suppressed by Landau 
damping. It, therefore, appears that a stationary reconnection process can operate in 
the turbulent current sheet only if cooling dominates heating in the reconnection region 
[20]. This is unlikely to be satisfied if the scale of the current sheet along the field lines is 
comparable with the size of the active region. This problem could, however, be avoided 
if the current sheet were disintegrated into numerous fragments in which the energy 
dissipation occurs. The field configuration in this case would, however, be essentially 
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FIGURE 1. Percolation of the electric current, Jo, through a network of resistors. Lq and ro are the 
inductivity and resistance of "normal" high conductive elements. R* and C* are the abnormal resistance 
and electrostatic double layers of turbulent low conductive elements 



different from that used in the traditional model, and the question about the reconnection 
mechanism in this situation remains open. 

It, therefore, turns out that both theoretical analysis and observations suggest that 
flares in solar atmosphere operate in a form of transition waves, which propagate through 
the system in a state of catastrophic equilibrium. The wave triggers the energy release 
process in local regions, and also provides a feed-back to the whole system. Next 
section shows that such a behavior is rather common and drives many puzzling processes 
observed on the Earth. 



FLARES AS CURRENT PERCOLATION THROUGH A RANDOM 

NETWORK OF RESISTORS 

Plasma instabilities and overheating of the turbulent current sheet lead to a disruption 
of the magnetic structure in active regions into a random network of "normal" and "tur- 
bulent" domains. A qualitative model of these domains can be constructed in terms of 
resistor elements (see Fig. 1). Under the condition of global current conservation (which 
is a natural consequence of high inductivity of the corona where the current sheets are 
located) the current is forced to percolate through this network of resistors [21]. Since 
conductivity of the resistors depends on the current density, j^, a random redistribution 
of the local currents in this network occurs. This redistribution is accompanied by turbu- 
lization of new elements as the current density in a local region increases over a critical 
value, jik > j cr . The system, therefore, switches into a dynamically equilibrium phase 
being governed by a non-linear feedback from one element to another through variations 
of the resistivity of the individual elements and currents redistribution. 



This system can be described in terms of the percolation in a random network of 
resistors in which the resistance of each element, R±, depends on the local current 
density j± as: 

' R± = Rq for j ik < j cr 

R lk = R t ~W 6 R for M>j cr , 1 ; 

where Rq is a very low "normal" resistance of "good" elements and R t is a very high 
"abnormal" resistance of "bad" elements in the network. Distribution of currents and 
voltage on the elements of the network is controlled by Kirchhoff laws. In particular, 
7k = Xjik = 7o due to the conservation of the total current, 7o, at any fc-section of the 
network while the voltage along any ^-circuit with the same initial and final points in the 
network is AV\ = E Jik^k- 

Note that percolation through random network of conductive elements (or fractal clus- 
ters of elements) has a wide range of applications [28, 3] from creation of nanomaterials 
with given conductive and optical properties to the study of spreading of forest fires, 
infection deceases and magneto spheric storms [32]. This process was studied in direct 
experiments (in superconductive ceramic samples [30] and conductive graphite paper 
with random holes [1 1]) in numerical simulations [10] and analytically [23]). However, 
a main tool to study percolating systems is computer simulation, taking into account the 
character of interconnection between the elements, dimension of the lattice and fractal 
dimension of clusters formed by conductive elements. 

The percolating systems have the following common properties: (i) Phase transition 
of the system at a certain value of critical parameter usually referred to as a percolation 
threshold manifests itself in drastic changes of the system properties (see e.g. [2]), 
(ii) Self-clustering of the network's elements into fractal-like patterns with a power- 
law spectrum of basic parameters, and (iii) Non-linear behavior of the system due to 
feedback from individual elements. 

One of the key properties of the network is the threshold-like character of the global 
resistance R neU which depends on the density of "bad" resistive elements, p\ = p\{Jq), 
and on the total current in the turbulent current sheet, as R net (7o — 7 cr )~ a . The flaring 
starts as R net — >■ <*>. This is expected if the density of the "bad" elements increases or/and 
To — > 7 cr , at which the phase transition of the system occurs. 

A particular property of the percolating random network is that the number of current 
clusters, N(x), depends on the network parameters (such as the size, volume, length, 
etc.) as N(x) x~ n . This indicates that the amplitude-frequency spectrum of solar flares 
within this approach must be a power-law, N(W) W~ m , and the same spectrum is 
expected for flares in red dwarfs of UV Cety-type. 

The same is valid for particles accelerated in clusters of the random network of resis- 
tors. The electric field generated during flaring events in the clusters can be evaluated as 

2 

E* = — = (10 2 - 10 3 )^^i, where a* = ^ (io 2 - 10 3 )c% is the conductivity 

a* ajoi 4^v ef f 

in the flaring region under condition j > j cr = nec Si . Here c S; is the ion sound speed and 
O>0e is the electron plasma frequency. The energy of particles accelerated by the electric 
field, £ ~ eEJ z , can be as high as 

£ ~ eE*L = (1 - 10) GeV x nf^ 2 (y^) , (2) 



where n% = n/10 8 cm~ 3 and Tj = T/10 7 K. Since the number of clusters of the length 
l z in the network region is N(l z ) « l~ k the expected energy spectrum of the accelerated 
particles is iV(£ = eEJ z ) oc £~ fc 5 where the value of k depends on the basic properties of 
the network. This indicates that the percolating network scenario can provide us with a 
natural explanation of the observed power-law spectrum of solar cosmic rays. 

CONCLUSIONS 

The pre-flaring fine structure of the magnetic field in active regions observed with space 
missions differs from the static "potential+force-free" field configuration assumed in 
the traditional scenarios of flares. It appears that the field configuration is in a dynami- 
cal rather than static equilibrium and resembles the configuration of the field with shear 
observed in TOKAMAKs. Plasma instabilities and overheating of the turbulent current 
sheet prevent the energy release process from operating steadily and disintegrate the cur- 
rent sheet into numerous turbulent and "normal" domains. The dynamical equilibrium of 
this system, which is controlled by the level of global current can be explained in terms 
of the current percolation through a random network of resistors. This makes it possi- 
ble to explain the power-law statistical properties of flares and micro flares and particle 
acceleration with the power-law energy spectrum. 
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